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We present closed analytic expressions for the order af triangle diagram 
contributions to the matrix elements of the singlet and non-singlet axial vector 
currents between the vacuum and a quark-antiquark state. We have calculated 
these vertex functions for arbitrary momentum transfer and for four different 
sets of internal and external quark masses. We show that both the singlet and 
non-singlet vertex functions satisfy the correct chiral Ward identities. Using 
the exact expressions for the finite axial vector form factors, we check the 
quality and the convergence of expansions at production threshold and for 
asymptotic energies. 
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1 Introduction 



This paper is the third of a series that deals with the computation of the electro- 
magnetic and neutral current form factors of heavy quarks Q at the two-loop level in 
QCD (UI21- Knowledge of these form factors puts forward the aim of a completely 
differential description of the electroproduction of heavy quarks, e + e~ — > 7*, Z* — > 
QQX, to order ot 2 s in the QCD coupling. This, in turn, will allow precise calculations 
or predictions of many observables, including forward-backward asymmetries, which 
are relevant for the (re) assessment of existing data on h quark production at the Z 
resonance or which will become important for the production of b quarks and t 
quarks at a future linear collider ||]. A detailed list of references can be found in 

In [TJ |2] we presented analytic results to order a 2 s of the heavy quark vector and 
axial vector vertex form factors for arbitrary momentum transfer, up to anomalous 
contributions. Here we close this gap and analyze the triangle diagram contributions 
to the axial vector vertex which exhibit the Adler-Bell- Jackiw ( AB J) anomaly E] . 
We use dimensional regularization and the prescription of |H] for implementing 
the Dirac matrix 75 in D 7^ 4, which is an adaptation of the method of j^]. In 
Section 2 we set up the notation and discuss chiral Ward identities which must be 
satisfied by the appropriately renormalized axial vector current of a heavy quark, 
and the flavor-singlet and non-singlet axial vector currents, respectively. In Section 
3 we give closed analytic expressions to order a 2 s for the anomalous (pure singlet) 
contributions to the axial vector form factors for four different sets of internal and 
external quark masses. First we present our results for spacelike squared momentum 
transfer s; then the form factors are analytically continued into the timelike region, 
and we determine also their threshold and asymptotic expansions. We check that 
the form factors satisfy the chiral Ward identities discussed in Section 2. In the 
context of this check we have calculated also the triangle diagram contributions to 
the pseudoscalar vertex function and the vertex function of the gluonic operator 
GG. In jTU] the triangle contributions to the bb vertex function of the non-singlet 
current i^^/y^t — bj^b were computed for massless b quarks, using the method 
of dispersion relations. Our chirality-conserving form factors for this specific mass 
combination confirm that result. In Section 4, we expand our exact expressions for 
the finite axial vector form factors at the production threshold and at asymptotic 
energies. Detailed comparison of this expansions with the full result allows us to 
assess the convergence properties of these widely used expansion techniques [TTJ H2] . 
We conclude in Section 5. 

2 Chiral Ward Identities 

We consider QCD with Nf — 1 massless and one massive quark. This includes the 
case of Nf = 6 with all quarks but the top quark taken to be massless. We analyze 
the following external currents: 

• The currents J^^x) = Q{x)j l _ l 'y5Q(x) and J5 = Q(x)j5Q(x), where Q denotes the 
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field associated with the massive quark. These currents are called axial vector and 
pseudoscalar current in what follows. 

• The flavor- singlet axial vector current J^(x) = Yl Q.{ x )l ij,15Q.{ x ) i where the sum 
extends over all quark flavors including Q. 

• The flavor non-singlet neutral axial vector current J^jf = 5^ =1 ^T/ZysTs^, where 
ipi = (u, d)i is the i-th generation quark isodoublet and r 3 is the 3rd Pauli matrix. 
This current couples to the Z boson. The corresponding pseudoscalar current is 
denoted by J§ s ■ 

The current J^jf is conserved in the limit of vanishing quark masses, while the 
currents Js M and are anomalous at the quantum level El • 

2.1 Axial Vector and Pseudoscalar Vertex Functions 

In the following we study the matrix elements of the axial vector and pseudoscalar 
vertex functions between the vacuum and an on-shell qq pair to order a 2 in the 
QCD coupling, where q denotes one of the Nf quark flavors. The bare 1-particle 
irreducible (1PI) vertex functions involving J 5/J and J5 are denoted by T g M and T q , 
respectively. The kinematics is J{p) — > q{j>\) + q(j>2), with p = p± + P2- 

T q (s)= l5 + A q (s) . (2) 

Following [3] it is useful to distinguish between two types of contributions to A g At 
and A q and, likewise, for the other vertex functions considered below: those where 
the axial vector (pseudoscalar) vertex 7^75 (75) is attached to the external quark line 
q are called type A. For the axial vector and pseudoscalar currents defined above 
these are non-zero only for q = Q. Type B contributions are called those where 
7^75, respectively 75, is attached to a quark loop. To order a 2 these are the triangle 
diagrams Fig. 1 (a) for the heavy quark currents Js iAt and J5, and the diagrams Figs. 
1 (a,b) for , , J$J, and J^ s . For massless quarks in the triangle Fig. 1 (b) 
only Fig. 1 (a) contributes to the matrix elements of J§ and J§ . 

The anomalous Ward identities, which will be discussed in Section 2.3, involve 
also the truncated matrix element of the gluonic operator GG between the vacuum 
and an on-shell qq pair. Here 

G(x)G(x) = e llva pG a ^ v {x)G a ' aP {x) , (3) 

with G a ^ u being the gluon field strength tensor. The operator © induces, to zeroth 
order in the QCD coupling, the momentum-space vertex Sk^k^e^ap with the gluon 
momentum k\ {k,2) flowing into (out of) the vertex. The truncated one-loop qq 
vertex function F q which will be needed in Section 2.3 is shown in Fig. 2. 

We use dimensional regularisation in D = 4 — 2e space-time dimensions. For 
the treatment of 75 we use a pragmatic approach: the calculation of the type A 
contributions is made using a 75 that anticommutes with 7^ in D dimensions. This 
approach works for open fermion lines, but is known to be algebraically inconsis- 
tent, in general, for closed fermion loops with an odd number of 75 matrices. The 
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Figure 1: Triangle diagram (type B) contributions to the various axial vector and 
pseudoscalar vertex functions considered in this paper. The double and single 
straight line triangles in (a) and (b) denote massive and massless quarks, respec- 
tively. The external q, q quarks are either massive or massless. The dashed line 
denotes either an axial vector or a pseudoscalar vertex. The crossed diagrams are 
not drawn. 
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Figure 2: One-loop anomalous vertex function F q . 



appropriate implementation of 75 for type B diagrams is the prescription of 't Hooft 
and Veltman Jj] or its adaption by [7J |H] where 7 

ry — l_ e ^Ml^M2^M3^4 (A) 

15 ,j c /iiAi2M3A 1 4 fill: V*J 

7 M 75 = 7^W™7 M W 3 (5) 

is used within dimensional regularisation. In these relations, the Lorentz indices 
are taken to be D-dimensional. The product of two e-tensors is thus expressed as a 
determinant over metric tensors in D dimensions. We employ these relations in this 
paper. 

In the pragmatic approach chiral invariance is respected, to second order in the 
QCD coupling, also by the unrenormalized type A contributions. For the above 
vertex functions the following well-known chiral Ward identity can be derived in 
straightforward fashion: 

Z^C = 2m ^ A Q ] - S Q(Pi)75 - 1^ Q (-P2) , (6) 

7 We use the convention e i23 = — e 0123 = +1. 
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where m 0i Q and Eg denote the bare mass and self-energy of Q. Because the deriva- 
tion of © was based on an anticommuting 75 we have 



- 7sSq(-p 2 ) = S Q (p 2 )7 5 . (7) 
2.2 Renormalization of the Currents 

We renormalize the ultraviolet (UV) divergences by defining a s in the MS scheme, 
while the renormalized quark mass and the quark wavefunction renormalization 
constant Z2 are defined in the on-shell scheme. That is, the renormalized quark 
self-energy is zero on-shell. For the on-shell mass tuq of the quark Q we have 
mo = Z° s mQ, where Z° s denotes the mass renormalization constant in the on-shell 
scheme. We described this hybrid renormalization procedure in detail in QH2]- 

Eq. (JBJ) implies that, within the pragmatic approach and to second order in 
the QCD coupling, neither the axial vector current J 5ai nor the operator ttiq^Q^Q 
do get renormalized with respect to type A contributions. Thus the renormalized 
second-order axial vector and pseudoscalar vertex functions of type A are given by 

r (A),R _ 7 OS v {A) ( o\ 
1 Q,n - Z 2 1 Q,» > l°J 

r (A),R _ os (A) _ 7 OS 7 OS v {A) / Q \ 

m Q L Q - m 0,Q Z 2 1 Q - m Q Z m Z 2 1 Q ■ { J > 

Nevertheless, the axial vector current must be renormalized because of the ABJ 
anomaly which is present in the triangle diagrams. When using an anticommuting 
75 for computing the type A and (jlj), (0) for the type B contributions the correctly 
renormalized current to order a 2 s turns out to be 

€ = ZjZj mUe J 5 , , (10) 

i.e., the second order triangle diagram contributions to the axial vector vertex func- 
tion are renormalized according to 

A qT = A qI + (*J + 5 5)7m75 , (11) 



where 

5 j EE Zj-l 



5 5 ee Z f j inite -1. (12) 



The renormalization constant Zj removes the UV divergences present in the dia- 
grams Figs. 1 (a,b). From their computation - see the next section - we find that 
this constant, which we define in the MS scheme, is given by 

Z J = l + 3 -^(^) 2 C\e) + 0(al), (13) 

where C F = (iV c 2 - 1)/2N C , T R = 1/2, C(e) = (4vr)T(l + e), and N c denotes the 
number of colors. This result is in agreement with jH]. The use of (HJ) and (0) 
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requires, in addition, a finite renormalization zj mtte in order to obtain the correct 
anomalous Ward identity for the axial vector current j^j- That is to say, for the 
renormalized (singlet) axial vector currents the non-renormalization J3j of the ABJ 
anomaly is to be maintained. For the axial vector current this reads 

(&>Jsp) R = 2im Q (Q l5 Q) R + ^-T R (GG) R . (14) 

From our calculation of the respective vertex functions involving massive and/or 
massless quarks, presented in Section 3, we find that Zj intte is given by 

7 finite _ , , 3C F T R (Us\ 2 n(n Z\ 



z J= 1+ -^r\^) ( 15 ) 

Again, we can compare with the analogous result obtained in [Sj for massless QCD, 
and we find agreement. 

For a pseudoscalar vertex the triangle diagrams Fig. 1 (a) are UV and infrared 
finite. This implies that, to order a 2 s , 

= Ajf > , (16) 

and mo&Q^Q = mQiQlzQ) r within the pragmatic approach. We emphasize that 
no finite renormalization of the pseudoscalar operator (Q75Q is necessary in order 
that the Ward identities Eqs. (J23|) and (J25J) below are satisified. In the existing 
literature, the issue of a potential finite renormalization of the pseudoscalar singlet 
operator was not resolved up to now jH]. Previous calculations [H] circumvented 
this issue by forming appropriate non-singlet combinations insensitive to this renor- 
malization. 

Next we consider the renormalization of the truncated QQ vertex function « s Fq(s) 
of the operator a s GG. This operator mixes with d^J^ under renormalization [15 : 

(GG) R = Z Gd GG + Z GJ , (17) 

where is the flavor-singlet current. The renormalization constants which we need 
here were given, e.g., in [S]. In the MS prescription we have 

Z Gd = 1 + 0(a s ) , (18) 
Z GJ = 6 -^^C(e) + 0(al). (19) 

Eqs. (fTTj) and (fT9"|) imply that the renormalized one-loop vertex function Fq is given 
by 

FX = F Q +iZ GJ p» ltil5 . (20) 

We recall that Eq. (jHJ is to be used for 7^75. 

Finally we consider the non-singlet currents J^f and J^ s , for definiteness in 
Nf = 6 flavor QCD with five massless quarks and the massive quark Q being the 
top quark. The divergence d^J^f = 2m 0i tt75t is non-anomalous, and within our 
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prescription for implementing 75 the operators J^f and mQ^^t do not get renor- 
malized. Let us consider the 1PI vertex functions involving J^jf and J^ s and an 
external on-shell ti pair. Using again the on-shell scheme for defining the wavefunc- 
tion and mass of the top quark, the respective renormalized vertex functions, which 
contain both type A and type B contributions, are 

rtr = z? a r%, (21) 

T NS,R = Z 0S Z 0S T NS_ (22) 

2.3 Renormalized Ward Identities 

Let us now discuss the chiral Ward identities for the various renormalized vertex 
functions introduced in Section 2.2. First to those of the axial vector and pseu- 
doscalar currents. For the renormalization scheme specified above one derives from 
Eq. @ that the type A contributions to the on-shell second order renormalized 
vertex functions (jSJ), © must satisfy the Ward identity 

p^ R = 2m Q T^ R , (23) 

In j2j we computed the vertex function Tq' R to order a 2 using an anticommuting 
75. We have checked by explicit computation of the pseudoscalar vertex function 
Tq^'^, also to second order in the QCD coupling, that our results given in [2] satisfy 
Eq. flUD. 

For the second order renormalized type B contributions the anomalous Ward 
identity 

P-A™ 8 = 2»0^' '-i*T H F§ (24) 

must hold. This equation follows from Eqs. (JUj) and (fTTj) . 

With Eq. ()24|) we can immediately write down the Ward identity for the triangle 
diagram contributions to the QQ vertex function of the flavor singlet current J§ . 
It reads 

N f 

E A qT( s > m i> m <?) = 2m Q A Q m ^ m «) - ^ N f T ^ F $( s > m <?) > ( 25 ) 

where rrij denotes the mass of the quark in the triangle loop. All triangle diagrams 
Figs. 1 (a,b) contribute to the left side of (J2HJ), while for the pseudoscalar vertex 
the Dirac trace over a triangle is non-zero only for the massive quark. 

Next we consider the vertex functions A^„ and A,? which correspond to the 
matrix elements of J^r and between the vacuum and an on-shell qq pair with 
q 7^ Q. To second order in the QCD coupling the only contributions that are non- 
zero are those of type B. For a massless on-shell qq pair we have A^ = to order ot 2 s 
due to the chiral invariance of the quark-gluon vertices. The same statement holds 
for the vertex function of (GG)r. Thus for external massless on-shell qq quarks 
the corresponding second order renormalized axial vector vertex function must fulfill 



P M A« = . (26) 
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The renormalization prescription is that of Eqs. (jSJ) and (jllj) . 

Finally to the non-singlet vertex functions. Because the operator equation 
d^J^f = 2mo,tt75t does not get renormalized within the pragmatic approach, the 
renormalized on-shell vertex functions ()21|) and (}2*2*|) must satisfy the canonical chiral 
Ward identity 

p » T xs,R = 2m t T? s ' R . (27) 

The type A and type B contributions to the vertex functions fulfill Eq. (j2*Tj) sep- 
arately. The type B contributions to r^f and T^ s are shown in Figs. 1 (a,b). 
Because we take massless u, d, c, s, and b quarks, only the t, b quark pair survives 
the cancellation of the weak isodoublets in the triangle loops. In fact Tf s = T t be- 
cause we have put m& = 0. If one takes a massless on-shell qq pair instead of ti, the 
corresponding non-singlet axial vector vertex function must vanish upon contraction 
with p^\ 

P^ R = . (28) 

In the next section we show by explicit computation that Eqs. (J23J, (jSEj) , (J2HJ) , (|2*7j) . 
and ()28)1 are satisified by the respective type B contributions. 



3 Results for the Triangle Diagrams 

We compute now the contributions Figs. 1 (a,b) to the axial vector and pseudoscalar 
vertex functions to order a 2 s and the one-loop anomalous vertex function F q . The 
on-shell vertex functions are sandwiched between Dirac spinors 

6(pi)Ag v(p 2 ) , (29) 
tiCpOAf) v(p 2 ) , (30) 
u( Pl )F q v(p 2 ), (31) 

where color indices are suppressed. 



3.1 Calculation 

Due to parity and CP invariance of QCD the following decompositions hold in four 
dimensions: 

A S = lfil5Gi(s,m i ,m q ) + ^—p tl j 5 G2{s,mi,rn q ), (32) 

Af) = A{s,m h m q ) 75j (33) 
F q = 2im q f(s,m q )'j 5 . (34) 

The form factors G±, G 2 , A, and / are dimensionless. As before the mass of the 
external quark is denoted by m q , while is the mass of the quark circulating in 
the triangle loop. 

Because of helicity conservation the form factors G 2 = A = F q = for massless 
external quarks. Moreover, A = when the quarks in the triangle of Fig. 1 (b) are 
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massless. Thus only if both the internal and external quarks are massive. 

We compute the renormalized type B two-loop vertex functions Figs. 1 (a,b) for the 
mass combinations (wij, m q ) = (m, m), (0, m), (m, 0), and (0, 0). This is adequate for 
a number of applications, for instance to Z* — > qq and Z* — > QQ. Furthermore, we 
determine f(s,m). Let us first consider the axial vector vertex. Its renormalization 
is specified in Eq. flll|) . The functions A q B f}' R are also infrared (IR) finite for D — > 4 
for any of the above mass combinations. Thus the renormalized form factors Gj tR 
are obtained from these functions by appropriate projections, which could be carried 
out, in principle, in four dimensions using ()32|). We consider the spinor traces [2] 

Vi, R = Tr [ 7m75 (fa + m q ) (Ag> + (5j + S 5 ) W ) (fa - m q )\ , (35) 

V 2 ,r = Tr [ 75 PlM (fa + m q ) (Ag + (5 j + 5 5 )y^ 5 ) (fa - m q )] . (36) 

Rather than calulating A q B ^' R first and then performing in (|33|) . (|3l)|) the index 
contractions and traces in four dimensions, it is of course technically much simpler 
to compute first the individual, UV-divergent pieces appearing on the right hand 
sides of (|33j) and (j3T?j) in D dimensions. Therefore 7^75 and 75 appearing in the 
above traces have to be treated according to (HJ), ©, even if two 75 appear in the 
same Dirac trace. Because "P^r, "P2,_r are finite quantities, the relations between 
them and the form factors can be evaluated in four dimensions. We obtain 

G 1 , R (s,m i ,m q ) = As{Aml _ s) , (37) 

and, if m q ^ 0, 

n ( v m q sV hR + (6m 2 - s)V 2 , R 

G 2jR (s,mi,m q ) = -m q g 2 (4m 2_ g ) ■ ( 38 ) 

The form factor associated with the triangle diagram contributions to the pseu- 
doscalar vertex function is both UV- and IR-finite and is given by 

1 r ( r>\ ~| 

A R (s,m,m) = -— Tr 75 (fa + m)lY Q ! (fa - m) . (39) 

Eq. (jlj is to be used in the computation of both the triangle loop Fig. 1 (a) and of 
the trace (J3HJ). 

The anomalous vertex function F q is nonzero only for a massive quark. The renor- 
malized form factor, which is IR-finite, is given by 

f R (s, m) = Tr [ 75 (fa + m) (Fq + iZ GJ p» llxl5 ) (fa - m)] . (40) 

The above form factors are computed with the technique that was applied in 
PUI2] to the two-loop vector and type A axial vector form factors. (Cf. these papers 
for detailed list of references.) Performing the D-dimensional 7 algebra using ^H] we 
obtain the form factors expressed in terms of a set of different scalar integrals. These 
integrals are reduced [T2j to a small set of master integrals. These integrals were 
evaluated with the method of differential equations [13111112011211 in [221I23I211- 
The master integrals and thus the above form factors are expressed in terms of 1- 
dimensional harmonic polylogarithms H(a; x) up to weight four [23 [2H] which are 
functions of the dimensionless variable x defined by Eq. (pfTj) below. 
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3.2 Results for Spacelike s < 

We give our results first in the kinematical region in which s is spacelike. Here the 
form factors are real. The form factors involving an internal and/or external mass 
to are also expressed in terms of the variable 

y/—S + Am? - v^-s 
x = = — 1—, (41) 
V — s + 4m 2 + v — s 

with < x < 1. Defining in four dimensions: 

G^{s,mi,m q ) = CFT R Qj(s,mi,m q ) (42) 

for j = 1,2 we obtain for the case (to, to): 

Gi(s,m,m) = 3 ln(TO 2 //i 2 ) + (32a; 2 ff(0, 0,1,0; x) -32a; 2 ff(l, 0,0,0; x) 

-64 a; 2 ff(0, 0, -1, 0; x) - 64 x 2 ff(0, 1, 0, 0; x) 
+4 (3 x 3 + x 2 - 4 C(2) x - 3 x - 1) x ff (0, 0; x) 

+4 (x 4 - 2 x 2 - 8 C(2) x 2 + 1) H(l, 0; x) )/ [(x - 1) (x + l) 3 ] 

+ ( 8 (2 x 3 - 3 x 2 + 4 x + 1) x 2 ^(0, 0, 0; x) 

+ (3x 6 -2x 5 + 8C(2)x 5 -8C(2)x 4 -3x 4 + 32C(3)x 4 

+ 32 C(2) x 3 - 64 C(3) x 3 + 4 x 3 + 32 C(3) x 2 - 8 C(2) x 2 - 3 x 2 

- 2 x + 8 C(2) x + 3) if (0; x) ) / [(x + l) 3 (x - l) 3 ] 

+ 16x (ff(0, 1,0; x) - 2 ff (1,0, 0;x) - 2 ff (0, -1, 0; x) ) / (x + l) 2 
-8 (x- l)ff(-l,0;x)/(x+ 1) 

+ 1/5 (-45 x 5 + 10 C(2) x 5 - 50 C(2) x 4 + 80 C(3) x 4 - 45 x 4 

- 80 C(3) x 3 + 90 x 3 - 40 C(2) x 3 + 16 C(2) 2 x 3 

- 16 C(2) 2 x 2 - 80 C(3) x 2 + 90 x 2 - 40 C(2) x 2 - 50 C(2) x 

+ 80 C(3) x - 45 x - 45 + 10 C(2)) / [(x - l) 2 (x + l) 3 ], (43) 

where /i is the mass scale associated with the MS renormalization of J 5j ^ and ((n) 
denotes the Riemann zeta function. The chirality-flipping form factors do not get 
renormalized. For the completely massive case we find: 

g 2 (s,m,m) = 32(2x 3 -x 2 + 8x + 3)x 3 ff(0,0,0;x)/[(x+l) 3 (x-l) 5 ] 
-256 x 3 ff (0, 0, -1, 0; x) / [(x - l) 3 (x + l) 3 ] 
+64x 2 (ff(0,l,0;x) - 2 ff (1,0, 0;x) 
-2 ff (0, -1, 0; x) )/ [(x - l) 2 (x + l) 2 ] 

-8 x 2 (ff (0, 0, 0, 0; x) +6 ff (0, 0, 0, 1; x) )/ [(x - l) 3 (x + 1)] 
+32, x 2 ( ff (1, 0, 0, 0; x) - ff (0, 0, 1, 0; x) )/ [(x - 1) (x + l) 3 ] 
+32 (x 2 - 6 x + 1) x 2 ff (0, 1, 0, 0; x) / [(x - l) 3 (x + l) 3 ] 
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+4 (3 x 6 + 8 C(2) x 5 - 12 C(3) x 5 - 2 x 5 + 8 C(2) i 4 -3x 4 
+ 32 C(3) x 4 - 40 C(3) + 4 x 3 + 64 C(2) x 3 -3i 2 
+ 32 C(3) x 2 + 8 C(2) x 2 - 12 C(3) x + 8 C(2) x - 2 x 
+ 3)x#(0;x)/[(x+ l) 3 (x- l) 5 ] 

+4 (13 x 5 - 13 x 4 - 4 C(2) x 4 - 30 x 3 - 20 C(2) x 3 + 20 C(2) x 2 
- 2 x 2 + 4 C(2) x + x - 1) x H(0, 0; x) / [(x - l) 4 (x + l) 3 ] 

+ 16 (x 2 + 2 C(2) x + 2 x + 1) x H(l, 0; x) / [(x - 1) (x + l) 3 ] 

+8x(#(0,l;x) -4#(-l,0;x))/[(x- 1) (x + 1)] 

+8/5 (l5 C(2) x 5 - 5 x 5 - 8 C(2) 2 x 4 - 55 C(2) x 4 + 40 C(3) x 4 

- 5 x 4 + 10 x 3 - 100 C(2) x 3 - 40 C(3) x 3 - 60 C(2) x 2 
+ 10 x 2 - 40 C(3) x 2 + 8 C(2) 2 x - 5 x - 35 C(2) x 

+ 40 C(3) x - 5 C(2) - 5) x / [(x - l) 4 (x + l) 3 ]. (44) 

For the case of massless internal quarks we obtain: 

</i(s, 0, m) = 3 ln(m 2 //i 2 ) + 8 x 2 H(0, 1, 0; x) / [(x + l) 3 (x - 1)] 
+4x 2 H(0, 0, 0; x) / [(x + l) 3 (x - 1)] 
+8x 2 #(0,0,l;x)/[(x+ l) 3 (x- 1)] 
+ 16 x 2 H(0, 1, 1; x) / [(x + l) 3 (x - 1)] 
-2 (x 2 + 3 x - 2) x H(0, 0; x) / (x + l) 3 
-4 (x 2 + 3 x - 2) x H(0, 1; x) / (x + l) 3 
-2(3x 2 -2x + 3)#(l;x)/(x+ l) 2 
-2(x 2 + l)#(l,0;x)/(x+ l) 2 

+ (-3 x 4 + 2 x 3 + 8 C(2) x 2 - 2 x + 3) H(0; x) / [(x + l) 3 (x - 1)] 
-4 (x 2 + l)#(l,l;x)/(x+ l) 2 

-(9 x 4 + 6 C(2) x 4 + 18 x 3 + 16 C(2) x 3 + 8 C(3) x 2 - 40 C(2) x 2 
- 18x+ 16C(2)x-9 + 2C(2))/[(x+ l) 3 (x- 1)], (45) 



g 2 {s,0,m) = 32x 3 #(0,l,0;x)/[(x- l) 3 (x + l) 3 ] 
+ 16 x 3 H(0, 0, 0; x) / [(x - l) 3 (x + l) 3 ] 
+32 x 3 H(0, 0, 1; x) / [(x - l) 3 (x + l) 3 ] 
+64x 3 #(0,l,l;x)/[(x- l) 3 (x + l) 3 ] 
-4 (x 3 + 5 x 2 - 3 x + 1) x H(0, 0; x) / [(x + l) 3 (x - l) 2 ] 
-8 (3 x 2 - 2 x + 3) x H(l; x) / [(x + l) 2 (x - l) 2 ] 
-8 (x 3 + 5 x 2 - 3 x + 1) x H(0, 1; x) / [(x + l) 3 (x - l) 2 ] 
+4 (-3 x 4 + 2 x 3 + 8 C(2) x 2 - 2 x + 3) x H(0; x) 
/[(x-l) 3 (x + l) 3 ] 
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-8 (x 2 + l)x#(l,0;x)/[(x + l) 2 (x - l) 2 ] 

-16(x 2 + 1) xH(l, l;x)/[(x + l) 2 (x- I) 2 ] 

-8 (x 4 + C(2) x 4 + 2 + 8 C(2) x 3 - 16 C(2) x 2 + 4 C(3) a; 2 

+ 8 C(2) x - 2 x - C(2) - 1) x / [(x - l) 3 (x + l) 3 ]. (46) 

For the case of massless external quarks we get: 

<?i(s, m,0) = 3 ln(-s//i 2 ) +2(5x 2 -8x + 7)x 2 #(0,0;x)/(x-l) 4 
-8x 2 H(0, 1,0; x) / {x- l) 4 
-4x 2 F(0,0,0;x)/(x- l) 4 
-2 (-3 x 3 + 4 C(2) x 2 + 9x 2 -9x-4 C(2) x 

+ 3 + 4C(2))x#(0;x)/(x-l) 4 
-16x#(l,0,0;x)/(x- l) 2 
-16x#(0, -1, 0; x) I (x - l) 2 
+2 (x + 1) (x 2 + 1) H(l, 0; x) / (x - l) 3 
+6#(l;x) -8(x+l)#(-l,0;x)/(x- 1) 
-(2 C(2) x 4 + 9 x 4 + 8 C(3) x 3 - 8 C(2) x 3 - 36 x 3 
+ 54 x 2 - 36 x + 8 C(2) x + 8 C(3) x + 9 
-2C(2))/(x-l)\ (47) 

and 

Q^s, 0,0) = 3 ln(- S/ y) -9 + 2C(2). (48) 

Moreover, we have 

g 2 ( s ,m,o) = g 2 (s, 0,0) = 0. (49) 

As already mentioned above the form factor associated with the pseudoscalar vertex 
function is non-zero only for the completely massive case. It is given by 

A R (s, m, m) = (|0 2 C F T R ( - 16 x 2 H(0, 0, 0; x) / [(x - l) 3 (x + 1)] 

+4 (x + C(2) x - C(2) + 1) xH(0, 0; x) / [(x + 1) (x - l) 2 ] 
+4 (3 C(3) x 2 - 4 C(2) x - 6 C(3) x + 3 C(3)) x H(0; x) 

/[(x-l) 3 (x+l)] 
+12x#(0,0,0,l;x)/[(x+l) (x- 1)] 
-8C(2)x//(l,0;x)/[(x+l)(x-l)] 
+2 x H(0, 0, 0, 0; x) / [(x + 1) (x - 1)] 
+8 x H(0, 0, 1, 0; x) / [(x + 1) (x - 1)] 
-8 x H(0, 1, 0, 0; x) / [(x + 1) (x - 1)] 
-8x#(l,0,0,0;x)/[(x+l) (x- 1)] 

+8/5 (2 C(2) x + 5 x - 2 C(2) + 5) C(2) x / [(x + 1) (x - l) 2 ]) . (50) 
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Finally the form factor associated with the anomalous vertex function Fig. 2 is: 

f R (s, m) = (g) 2 C F ( 3 ln(m 2 //i 2 ) - (x - 1) ff(0, 0; x) / (x + 1) 
-2(x- l)#(0,l;x)/(x + 1) - (7x + 4C(2)x + 7-4C(2))/(x+ 1)). (51) 

3.3 Analytical Continuation above Threshold 

We perform the analytical continuation of the form factors to the physical region 
above threshold, s > 4m 2 , — 1 < x < 0, respectively s > for massless external 
quarks, with the usual prescription s — > s + iS. For the cases involving masses the 
variable x becomes a phase factor if < s < 4m 2 [I]. For s > 4m 2 we define: 

y= ^-^~^\ (52) 
y/s + v s — Am 2 

and the continuation in x is performed by the replacement: 

x — ► — y + id . (53) 

The real and imaginary parts of the form factors are defined through the relations: 

G jjR (s + iS) = Re G jtR (s) + m Im G j<R (s) , (54) 

and likewise for A R and f R , where also a factor 7r is taken out of the respective 
imaginary part. 

In the completely massive case we get for s > 4m 2 : 

Re&(s,m,m) = -8 (2y 3 + 3y 2 + Ay - 1) y 2 H(0, 0, 0; y) / [(y - l) 3 (y + l) 3 ] 
+32 y 2 H(-l, 0, 0, 0; y) / [{y + 1) (y - l) 3 ] 
+64 y 2 H(0, -1, 0, 0; y) / [(y + l)(y- l) 3 ] 
-192 C(2) y 2 H(0, -1; y) / [{y + 1) (y - l) 3 ] 
-32 y 2 H(0, 0, -1, 0; y) / [(y + 1) (y - l) 3 ] 
+64y 2 i/(0,0,l,0;i/)/[( Z /+l)( Z /-l) 3 ] 
-32 y #(0,1,0; 2/)/ (2/ -l) 2 
+16j/i7(0,-l,0; 2 /)/(j/-l) 2 
+4 (3 j/ 3 - j/ 2 - 3 j/ - 4 C(2) J/ + 1) y H(0, 0; y) 

/[(i/+l)(i/-l) 3 ] 
-32 y #(-1,0,0; y)/ (y - l) 2 
+96C(2)2/#(-l;2/)/(2/-l) 2 

+ (3 + 2 + 40 C(2) y 5 + 32 C(3) y 4 + 64 C(2) y 4 - 3 y A 
+ 64 C(2) y 3 + 64 C(3) y 3 - 4 y 3 + 32 C(3) y 2 - 32 C(2) y 2 
- 3 y 2 + 2 y - 8 C(2) y + 3) #(0; y) / [(y - l) 3 (y + l) 3 ] 

+8 (y + l)F(l,0;y)/(y-l) 
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-4 (y 4 - 2 y 2 + 16 C(2) y 2 + 1) H(-l, 0; y) / [(y + l)(y- if] 
-1/5 (45 y 5 + 170 C(2) y 5 + 80 C(3) y 4 + 70 C(2) y 4 - 45 y 4 

- 90 y 3 - 256 C(2) 2 y 3 - 200 C(2) y 3 + 80 C(3) y 3 

- 160 C(2) y 2 - 256 C(2) 2 y 2 + 90 y 2 - 80 C(3) y 2 

- 80 C(3) y + 45 y + 110 C(2) y - 45 + 10 C(2)) 
/[(y+lf(y-lf] 

+3 ln(m 2 //i 2 ), (55) 



lmg 1 {s,m,m) = -8 (2y 3 + 3y 2 + 4y - 1) y 2 #(0, 0; y) / [(y - l) 3 (y + l) 3 ] 
+32y 2 tf(-l,0,0;y)/[(y+l) (y - l) 3 ] 
+64y 2 ff(0,-l,0;y)/[(y+l) (y - l) 3 ] 
-32y 2 H(0,0,-l;y)/[(y+l)(y-lf] 
+Q4y 2 H(0,0,l;y)/[(y+l)(y-lf] 

+4 (3 y 3 - y 2 - 4 C(2) y - 3 y + 1) y if (0; y) / [(y + 1) (y - if] 

-32yff(-l,0;y)/(y-l) 2 

+ 16yH(0,-l;y)/(y-lf 

-32ytf(0,l;y)/(y-l) 2 

-4(y + l)tf(-l;y)/(y-l) 

+8 (y + l)tf(l;y)/(y-l) 

+ (3y 4 -4y 3 + 8C(2)y 3 + 2y 2 + 32C(3) y 2 -4y-8((2)y 
+ 3)/[(y + l)(y-l) 3 ]. (56) 

For the chirality-flipping form factor we have 

Re£ 2 (s,m,m) = 32 (2 y 3 + y 2 + 8 y - 3) y 3 H(0, 0, 0; y) / [(y - l) 3 (y + l) 5 ] 
-256y 3 tf(0,0,l,0;y)/[(y+l) 3 (y-l) 3 ] 
-384C(2)y 2 tf(-l;y)/[(y+l) 2 (y-l) 2 ] 
-8y 2 tf(0,0,0,0;y)/[(y+l) 3 (y-l)] 
-32 (y 2 + 6 y + 1) y 2 H(0, -1, 0, 0; y) / [(y + l) 3 (y - l) 3 ] 
+ 128 y 2 H(-l, 0, 0; y) / [(y + l) 2 (y - l) 2 ] 
-32 y 2 H(-l, 0, 0, 0; y) / [(y + 1) (y - if] 
-64y 2 #(0,-l,0;y)/[(y+l) 2 (y-l) 2 ] 
+ 128y 2 tf(0,l,0;y)/[(y+l) 2 (y-l) 2 ] 
+32y 2 if(0,0,-l,0;y)/[(y+l)(y-l) 3 ] 
+48 y 2 H(0, 0, 0, -1; y) / [(y + l) 3 (y - 1)] 
-4 (13 y 5 - 2 C(2) y 4 + 13 y 4 - 30 y 3 - 14 C(2) y 3 - 14 C(2) y 2 
+ 2 y 2 - 2 C(2) y + y + 1) y H(0, 0; y) / [(y + l) 4 (y - l) 3 ] 
+8 (y 4 + 12 C(2) y 3 - 2 y 2 + 72 C(2) y 2 + 12 C(2) y 
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+ l)yiJ(0,-l;y)/[(y+l) 3 (y-l) 3 ] 
+ 16 (y 2 + 4 C(2) y - 2 y + 1) y H(-l, 0; y) / [(y + 1) (y - l) 3 ] 
-4 (3 y e + 12 C(3) y 5 + 40 C(2) y 5 + 2 y 5 + 32 C(2) y 4 

+ 32 C(3) y 4 - 3 y 4 + 40 C(3) y 3 - 4 y 3 + 128 C(2) y 3 
+ 32 C(3) y 2 - 3 y 2 - 64 C(2) y 2 + 2 y + 12 C(3) y 
- 8 C(2) y + 3) y if (0; y) / [(y - l) 3 (y + l) 5 ] 
-32yff(l,0;y)/[(y+l)(y-l)] 
+4/5 (10 y 4 + 165 C(2) y 4 + 80 C(3) y 3 + 29 C(2) 2 y 3 

- 80 C(2) y 3 - 20 y 3 - 314 C(2) 2 y 2 - 170 C(2) y 2 

+ 29 C(2) 2 y + 80 C(2) y + 20 y - 80 C(3) y - 10 

+ 5C(2))y/[(y + l) 3 (y-l) 3 ], (57) 



and 



lmg 2 (s,m,m) = -256 y 3 ff(0, 0, 1; y) / [(y + l) 3 (y - l) 3 ] 

+32 (2 y 3 + y 2 + 8 y - 3) y 3 ff(0, 0; y) / [(y - l) 3 (y + l) 5 ] 

+128y 2 ff(0,l;y)/[(y+l) 2 (y-l) 2 ] 

+32y 2 F(0,0,-l;y)/[(y+l)(y-l) 3 ] 

-64y 2 ff(0,-l;y)/[(y+l) 2 (y-l) 2 ] 

-32 (y 2 + 6 y + 1) y 2 ff(0, -1, 0; y) / [(y + l) 3 (y - l) 3 ] 

+ 128y 2 /f(-l,0;y)/[(y+l) 2 (y-l) 2 ] 

-8y 2 tf(0,0,0;y)/[(y+l) 3 (y-l)] 

-32y 2 //(-l,0,0;y)/[(y+l) (y - l) 3 ] 

-32yff(l;y)/[(y+l)(y-l)] 

+16yff(-l;y)/[(y + l)(y-l)] 

-4 (13 y 5 + 13 y 4 + 2 C(2) y 4 - 30 y 3 - 18 C(2) y 3 + 2 y 2 

-18C(2)y 2 + y + 2C(2)y+l)yiJ(0;y) 

/[(y+l) 4 (y-l) 3 ] 
-4 (3 y 4 + 12 C(3) y 3 - 4 y 3 + 8 C(2) y 3 + 2 y 2 + 8 C(3) y 2 

- 4 y + 12 C(3) y - 8 C(2) y + 3) y / [(y + l) 3 (y - l) 3 ] . 

(58) 

We obtain for the non-zero pseudoscalar and anomaly form factors: 

ReA R (s,m,m) = (^f S^f^yjs [4 (3 C(3) + 3 C(3)y 2 + 6yC(3) - 8yC(2))ff (0; y) 

-2 (1 + y) (2 - 2 y + y((2) + C(2))ff(0, 0; y) + 16 yH(0, 0, 0; y) 
-12 (1 + y) 2 /f (0, 0, 0, -1; y) + 2 (1 + y) 2 /f (0, 0, 0, 0; y) 
-8 (1 + y) 2 /f (0, 0, -1, 0; y) - 24((2)(1 + y) 2 /f (0, -1; y) 
+8 (1 + y) 2 ff (0, -1, 0, 0; y) - 16 C(2)(l + y) 2 H(-l, 0; y) 



14 



+8 (1 + y) 2 H(-l, 0, 0, 0; y) - ^ y( 2 (2) - 4 y 2 C(2) 
-f C 2 (2)+4C(2)-|C 2 (% 2 



lmA R (s,m,m) = — 



2tt7 (i- y )(i + y y 



2(l + y)(yC(2) + 2y-2 + C(2))//(0;y) 



+ 16 y#(0, 0; y) + 2 (1 + y) 2 #(0, 0, 0; y) - 8 (1 + y) 2 #(0, 0, -1; y) 
+8(l + y) 2 #(0,-l,0;y) 

+8 (1 + y) 2 H(-l, 0, 0; y) + 24 y((3) + 12 C(3)y 2 + 12 C(3)l , (60) 



Ref R (s,m) 



2C f 



,2tt/ 1-y 

+3 (1 - y) ln(m 2 //. 2 ) - 7(1 - y)] , 
lmf R (s,m) = ^^L(l + y )H(0;y). 



(1 + y) if (0, 0; y) - 2(1 + y) if (0, -1; y) + C(2)(l + y) 

(61) 



In the case of massless internal quarks we get for s > Am 2 : 
ReC/i(s,0, m) 



(62) 



-8y 2 tf(0,0,-l;y)/[(y-l) 3 (y + l)] 
+4y 2 /7(0,0,0;y)/[(y-l) 3 (y+l)] 
-8y 2 ff(0,-l,0;y)/[(y-l) 3 (y + l)] 
+ 16y 2 ff (0, -1, -l;y) / [(y - i f (y + 1)] 
+4(y 2 -3y-2)yff(0,-l;y)/(y-l) 3 
-2 (y 2 -3y-2)yff(0,0;y)/(y-l) 3 
-(3y 4 + 2y 3 + 4((2) y 2 - 2y - 3) H(0;y) 

/[(y-l) 3 (y + l)] 
+2 (3y 2 + 2y + 3)ff(-l;y)/(y-l) 2 
-4(y 2 + l)ff(-l,-l;y)/(y-l) 2 
+2 (y 2 + l)ff(-l,0;y)/(y-l) 2 

-(9 y 4 - 18 y 3 - 4 C(2) y 3 - 10 C(2) y 2 + 8 C(3) y 2 + 18 y 

-4C(2)y + 2C(2)-9)/[(y-l) 3 (y + l)] 
+3 ln(m 2 //i 2 ), (63) 



Im^( S ,0,m) = -8y 2 fJ(0,-l;y)/[(y-l) 3 (y + l)] 
+4y 2 ff(0,0;y)/[(y-l) 3 (y+l)] 
-2 ( y 2 -3y-2)yff(0;y)/(y-l) 3 
+2(y 2 + l)ff(-l;y)/(y-l) 2 

-(3 y 4 + 2 y 3 - 4 C(2) y 2 - 2 y - 3) / [(y - l) 3 (y + 1)], (64) 
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Re£ 2 (s,0,m) = 32 y 3 H(0, -1, 0; y) / [(y + l) 3 (y - l) 3 ] 



+32y 3 tf(0,0,-l;y)/[(y+l) 3 (y-l) 3 ] 

-64 y 3 #(0, -1, -1; y) / [(y + if (y - l) 3 ] 

-16y 3 #(0,0,0;y)/[(y+l) 3 (y-l) 3 ] 

+4 (3 y 4 + 2 y 3 + 4 C(2) y 2 - 2 y - 3) y H(0; y) 

/[(y + l) 3 (y-l) 3 ] 
-8 (3y 2 + 2y + 3) y F(-l; y) / [(y - l) 2 (y + l) 2 ] 
-8 (y 3 - 5 y 2 - 3 y - 1) y H (0, -1; y) / [(y - l) 3 (y + l) 2 ] 
-8 (y 2 + 1) y (-1, 0; y) / [(y - l) 2 (y + l) 2 ] 
+4 (y 3 - 5 y 2 - 3 y - 1) y H(0, 0; y) / [(y - l) 3 (y + l) 2 ] 
+16 (y 2 + 1)2/^-1, -1; y) / [(y - l) 2 (y + l) 2 ] 
-4 (C(2) y 4 - 2 y 4 + 4 y 3 + 4 C(2) y 3 + 8 C(2) y 2 - 8 C(3) y 2 

- 4y + 4C(2) y + 2 - C(2)) yj [(y + l) 3 (y - l) 3 ], (65) 



Im0 2 (s,O,m) = 32y 3 #(0,-l;y)/[(y+l) 3 (y-l) 3 ] 



-16y 3 tf(0,0;y)/[(y+l) 3 (y-l) 3 ] 

+4 (y 3 - 5 y 2 - 3 y - 1) y H(0; y) / [(y - l) 3 (y + l) 2 ] 

-8 (y 2 + l)ytf(-l;y)/[(y-l) 2 (y+l) 2 ] 

+4 (3 y 4 + 2 y 3 - 4 C(2) y 2 - 2 y - 3) y / [(y + l) 3 (y - l) 3 ]. 



Now to the case of massive internal and massless external quarks. An on-shell two- 
gluon intermediate state in Fig. 1 (a) could lead for s > to an imaginary part of 
G 2 only (for an off-shell Z boson), as a consequence of the Landau- Yang theorem 
[2*7] . But this form factor is zero for massless on-shell q, q. The form factor G\ has 
an absorptive part only if s > 4m 2 . In the region < s < Am 2 G\ is real and we 
get: 



denotes the Clausen function of second (third) order [2*5] . 
For s > Am 2 we find: 



Re^( S ,m,0) = 3 1n( S //i 2 ) + 2 (5y 2 + 8y + 7) y 2 #(0, 0; y) / (y + l) 4 



and 



(66) 




4 166 6 166 

+—Cl 2 (26) + -^-Cl 2 {26) + -Cl 2 {26) + —^Cl 2 (A6) 
rw r z w r 

--Cl 2 (A6) - -Ck(26) + ^Ch(26) + -C7 3 (40) , (67) 



where r = s/m 2 , w 
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and 



-4y 2 H(0,0,0;y)/(y+l) 4 
+8y 2 H(0,-l,0;y)/(y+l) 4 
+2(3y 3 + 9 2 / 2 + 4C(2) 2 / 2 + 10C(2) 2 / + 9y + 4C(2) 

+ 3)yH(0;y)/(y+l) 4 
-16yH(-l,0,0;y)/(y + l) 2 
-16yH(0,l,0 ]y )/(y+l) 2 
-6 (y 2 - 8 C(2) y + 2 y + 1) H(-l; y) / (y + l) 2 
+8 (y-l)#(l,0;jO/(y+l) 
-2(y-l)(y 2 + l)tf(-l,0;y)/(y+l) 3 
-(9 + 32 C(2) y 4 + 36 y 3 + 56 C(2) y 3 - 8 C(3) y 3 
+ 42 C(2) y 2 + 54 y 2 + 36 y - 8 C(2) y - 8 C(3) y 
+ 9-2C(2))/( Z/ +l) 4 , (68) 

bn&(s,m,0) = 8j/ 2 #(0,-1;j/)/(j/ + 1) 4 
-Ay 2 H(0,0;y)/(y+l) 4 
+2 (5j/ 2 + 8y + 7) y 2 H(G;y)/{y + iy 
-16yH(0,l;y)/(y + lf 
-16yH(-l,0;y)/(y+l) 2 
+8 (y-l^l^/fo + l) 
-2 (y-lH^ + l^-l^/^ + l) 3 
+2 (3 2 / 3 + 9y 2 + 4C(2) 2/ 2 + 92/ + 6C(2)2/ + 3 

+ 4C(2))y/( 2 / + l) 4 -3. (69) 

For the completely massless case we have for s > 0: 

Re & (s, 0,0) = 3 ln(s//i 2 ) -9 + 2C(2), (70) 

and 

Im0i(s,O,O) = -3. (71) 

3.4 Check of the Chiral Ward Identities 

The anomalous Ward identity ()24|) reads in terms of the form factors: 



S Oi 

Gi tR (s, m, m) + - — -G 2 .r{s, m, m) = A R (s, m, m) + —^-T R f R (s, m) , (72) 
Am z 47r 

and Eq. (|25J) yields 

s ot 
G 1)R (s, 0, m) + ^-^G 2 , R (s, 0, m) = -^T R f R (s, m) . (73) 

The results of Section 13.31 satisfy this equation, both for spacelike and timelike s. 
Eqs. (j^Sj) and (J25j) are trivially satisfied because G\ iR u{p\) ^75^(^2) = for massless 
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external quarks. 

For the type B contributions the non-singlet Ward identity (|27|) has the form: 

AG 1>R + -^AG 2ifi = A R {s, m, m) , (74) 

where 

AGj^R = Gj t n(s, m, m) — Gj t ji(s, 0, m) = Gj(s, m, m) — Gj(s, 0, m) . (75) 

As Eqs. f!72j) and ([7H|) are satisfied by our form factors, Eq. (|73j) is of course fulfilled, 
too. 

A specific combination of the above vertex functions, namely Gi(s,m, 0) — 
Gi(s, 0,0), which is relevant for Z — > massless bb, was computed in [TU] using 
Cutkosky rules and dispersion relations. Applying the notation of that paper, one 
obtains the relation Gi(s, m, 0) — Qi(s, 0, 0) = —I2 where the function I 2 is given in 
[TU] for s < and < s < Am 2 . Using in the respective kinematical regions the real 
and imaginary parts of the form factors Q\ given above, an analytical comparison 
with Eqs. (2.27), (2.28), (2.30) and (2.31) of [JO] yields full agreement. 



4 Threshold and Asymptotic Expansions 



Next we expand the axial vector and pseudoscalar form factors for a massive external 
QQ pair near threshold (s ~ 4m 2 ), in powers of f3 = yl — ^7^- This is useful for 
applications to, e.g., e + e~ — ► tt near threshold. 

Up to and including terms of order (j3) 2 - respectively of order (fl) 3 if the leading 
term vanishes - the real and imaginary parts are: 



Re Qi(s, m, m) 



Im Qi(s, m, m) 
Re G-zis, m, m) 



Im Q2(s, m, m) 
Re G\(s, 0, m) 

Im Gi(s,0 : m) 



1 Q 23 

-- C(2) - - + 16 C(2) ln(2) + 3 ln(m 2 //i 2 ) 



+(3 2 
P 3 



g C (2)-fc(2)ln(2) + |' 



32 , /n . 16 , . . 44 
- \n(P) + T ln(2) - — 



21 34 2 

— C(3) - 4 C(2) ln(2) + 6 Y C(2) + - + 4 ln(2) - 12 C(2) (3 



226 35 4 44 

- 6 - ^ C(2) + f C(3) + 3 ln(2) + - C(2) ln(2) 



-2 + 3C(2) + 4/3 



ln(2) 



-3-5C(2) 



4 8 23 

- C(2) - - ln 2 (2) + 8 ln(2) - - + 3 ln(m 2 //i 2 ) 



15 



ln 2 (2) 



38 4 



-4+3 ln(2) 
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15 15 
(3 2 ln(2) 



C(2) 
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ReQ 2 (s, 0, m) 



-^C(2) + ^ln 2 (2)-4 1n(2) + ^ 



+P 2 



38 16 , 16 , , ns 8 

-i5-y" l(2) + T ln(2) + 5 c(2) 



Im 6 2 (s, 0, m) = 2 - - ln(2) + f? 



-§+T-0> 



The form factor Qi(s, m, 0) is given, for r = s/m 2 — > 0, approximately by 



m, 0) = 3 ln(m ) + - - -r - — r 



while for s ~ 4m 2 its imaginary part has the small /3 expansion 
lm£i(s,m,0) = 0(s-4m 2 )/5 3 



16 64 
T HP) + 8 ln(2) - — 



(76) 



(77) 



(78) 



For brevity we give for the pseudoscalar form factor only the terms of order (/3)°: 



ReA R (s, m, m) = 

Im4(«,m,™) = (^) 3((2)C F T R . 



2 



5C(2)-yC(3) + 12C(2)ln(2) 



(79) 



Finally we study the behaviour of the form factors involving the massive quark 
for large squared momentum transfer s ^> m 2 . We obtain, up to and including 
terms of order r~ 2 , with L = ln(r): 



Re Qi(s, m, m) 



3 ln( S/ y) -9 + 2C(2) 

+ - (-L 2 + (3 - 4 C(2)) L — I + 12 C(2) - 8 C(3)) 
r 

2 + 24 L + 140 C(2) + 32 L((3) - 32 Z<(2) - 18 L 2 



+ 



4L 3 



256 



C 2 (2) + 8 C(2)L 2 + 8 (—2 L + 2)C(2) - 64 C(3) 



Im&(s,m,m) = -3 + - (2 L - 3 + 4 C(2)) 



1 



24 + 36 L - 4 L 2 - 32 C(3) + 32 C(2) - 16 L((2) 



Re (/2(s, m>, m ) 



Im 02(s, m, m) 
Re Qi(s, 0, m) 



-(L 2 -6L + 4-2C(2)) 

+ 4 — 32 L + 4 L 2 — 72 C(2) - 48 L((3) + 64 C(3) 

-H^C 2 (2) + ^ 4 -4C(2)L 2 + 32LC(2) 
- (3 - L) + \ [32 - 8 L + 48 C(3) - 8 L((2) — - L 3 — 32 C(2) 



3 \n(s/fi 2 ) - 9 + 2 C(2) + - (-L 2 + 3 L - 1) 
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Im Gi(s, 0, m) 
Re^ 2 (-s, 0, m) 
lmQ 2 (s, 0,m) 
Re Qi(s, m, 0) 



Im Gi(s, m, 0) 



+- 



- + 33 L - 10 C(2) + - L 3 + 8 C(3) - 4 L£(2) - 13 1} 
2 3 

2 1 r 

33 + 26L-2L 2 -4C(2) 



-3 + -(2L-3) + - 



r 



-(L 2 -6L + 4-2C(2)) + ^(-32L + 8C(2) + 12L 2 ) ; 
^(3-L) + -^(32 -24L), 



-9 + 3 ln( S /^ 2 ) + 2 C(2) + - C(3) - C(2) L 

r 



+- 



-3 



f -6C(2) + ^L 3 + 3L-L 2 



+ 12LC(2) + 8(2-2L)C(2) 
, »C(2) , 1 



16C(3) 



-3 + 2L-4C(2) -2L 2 ). 



(80) 



The leading asymptotic terms of the pseudoscalar form factor are 



Re A R (s, m, m) 



Im Ar(s, m, m) 



a 
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^ ) 2 C F T R l - [1/12 L 4 + 4 C(2) - y C 2 (2) 
2L 2 -C(2)L 2 -12LC(3)], 

ry \ 2 j 

-i) C F T fi -[4L + 12C(3)-2LC(2)-l/3L 3 ] 

Z7T/ r 



(81) 



The limit r — > oo corresponds to the massless limit m — > 0. Therefore the chirality- 
flipping form factors Q2 and expected, of order 1/r, and the terms of 

order r° in the above expansion of the chirality-conserving form factors Q\ are equal 
to Re(/i(s, 0, 0) and Im£?i(s, 0, 0), respectively. 

In Figs. 3 - 6 we have plotted the real and imaginary parts of the form factors 
CfTrQj as functions of the rescaled center-of-mass energy y / s/(2m), for \i — m 
and N c = 3. Also shown are the values of the asymptotic expansions of these form 
factors given in Eq. (|%U|) . and of the threshold expansions, Eqs. (I77j) and lj7B]). As 
far as the threshold expansions are concerned, these Figures show that both the 
next-to-leading terms in the expansions and restriction to a relatively small region 
above threshold are required in order to get a satisfactory approximation to the exact 
result. An exception is the expansion of Re£?i(s, m, 0) around s ~ where the terms 
up to second order approximate this form factor very well up to ^Js ~ 1.5m. The 
asymptotic expansions (JHUj). on the other hand, provide a very good approximation 
to the form factors over a wide range of center-of-mass energies. 
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Figure 3: The real and imaginary parts of CFT R Q 12 (s,m,m) as functions of the 
rescaled cm. energy and their threshold and asymptotic expansions. 
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Figure 4: The real and imaginary parts of CfTrQ\^{s^ 0, m) as functions of the 
rescaled cm. energy and their threshold and asymptotic expansions. 
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Figure 5: The real and imaginary parts of CfT R Qi(s,Tn,0) as functions of the 
rescaled cm. energy and their threshold and asymptotic expansions. 
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Figure 6: The real and imaginary parts of CfTrQi(s, 0, 0) as functions of the rescaled 
cm. energy. 



5 Summary 

In this paper we calculated the triangle diagram contributions to the form factors of 
a heavy-quark axial vector current, of the flavor- singlet and non-singlet axial vector 
currents, and of the corresponding pseudoscalar currents to second order in the QCD 
coupling, for four different combinations of internal and external masses. Moreover, 
we have given the threshold and asymptotic expansions of these form factors to a 
degree of approximation which should prove useful for applications. We compared 
our exact results with these expansions at threshold and for large energies. For the 
threshold expansions we found that they approximate the exact results only over 
a very restricted energy range above threshold, and that at least the second order 
terms in the expansion are required for a reliable description. The large energy 
expansions appear to converge very well, and do in general yield a good description 
even for energies well below the asymptotic regime. 

In this calculation we have used the implementation of the 75 matrix according 
to |7J E] in dimensional regularization. This method is algebraically consistent and 
convenient for computations, but breaks chiral invariance in D 7^ 4. To restore chiral 
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invariance in this scheme, finite counterterms are required. With the exception of 
the counterterm for the pseudoscalar singlet current, these were determined in jS] for 
massless quarks. In this work, we calculated all finite counterterms required by chiral 
invariance in this scheme for the case for massive quarks, including the counterterm 
for the pseudoscalar singlet current (which turns out to be zero). All counterterms 
agree with the massless results of [S], illustrating their purely ultraviolet origin. 
Using these counterterms, we have shown that the appropriately renormalized vertex 
functions satisfy the correct chiral Ward identities. This exemplifies for a non-trivial 
two-loop case that the counterterms of [H] can be applied, as expected, also to 
infrared-finite Green functions involving massive quarks. However, it is a separate 
issue to apply the method of jHj to Green functions (with massless and/or massive 
quarks) that are infrared-divergent - using dimensional regularization as well - and 
ask whether these Green functions satisfy the correct chiral Ward identities. We 
shall report on this issue and on applications of our results in this paper and in [TJ |2| 
to the electroproduction and, in particular, to the forward-backward asymmetry of 
heavy quarks in a future publication. 
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